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ABSTRACT 

We investigate pure luminosity evolution models for early-type (elliptical and 
SO) galaxies (i.e., no number density change or morphology transition), and exam- 
ine whether these models are consistent with observed number counts in the B, I and 
K bands and redshift distributions of two samples of faint galaxies selected in the I 
and K bands. The models are characterized by the star formation time scale tsf and 
the time t gw when galactic wind blows in addition to several other conventional pa- 
rameters. We find the single-burst model (tsf = 0.1 Gyr and t gw — 0.353 Gyr), which 
is known to reproduce the photometric properties of early-type galaxies in clusters, is 
inconsistent with redshift distributions of early-type galaxies in the field environment 
due to overpredictions of galaxies at z > 1.4 even with strong extinction which is at 
work until t gw . In order for dust extinction to be more effective, we change tsf and t gw 
as free parameters, and find that models with tsf ^ 0.5 Gyr and t gw > 1.0 Gyr can be 
made consistent with both the observed redshift distributions and number counts, if 
we introduce strong extinction (E(B — V) > 1 as a peak value). These results suggest 
that early-type galaxies in the field environment do not have the same evolutionary 
history as described by the single-burst model. 

Key words: cosmology: observations — galaxies: elliptical and lenticular, cD — 
galaxies: evolution — galaxies: formation — galaxies: photometry — galaxies: statis- 
tics 



1 INTRODUCTION 

How early-type galaxies were formed and evolved is a key is- 
sue in extragalactic astronomy which remains controversial. 
Evolution of early-type galaxies in clusters seems to be well 
expressed by the so-called single-burst model in which galax- 
ies experience a star burst at the initial phase of their forma- 
tion and then evolve passively without any subsequent star 
formation. For example, Ellis et al. (1997) found that the 
scatter in the colour-magnitude relation of early-type galax- 
ies is quite small for z ~ 0.5 clusters, indicating that the ma- 
jor star formation ended at z Z 3 (See also Stanford, Eisen- 
hardt & Dickinson 1998). Kodama et al. (1998) investigated 
the slope and zero-point of the colour-magnitude relation in 
17 distant clusters with redshift 0.31 < z < 1.27, and con- 



strained the epoch of major star formation to z > 2 — 4. How- 
ever, it is fairly controversial whether the single-burst model 
holds for early- type galaxies in the field environment. Abra- 
ham et al. (1999) analysed the spatially resolved colours of 
galaxies detected in the Hubble Deep Field (HDF; Williams 
et al. 1996) and argued that 40 % of early-type galaxies 
have experienced the major star formation at z < 1. Ko- 
dama, Bower & Bell (1999) obtained similar results on the 
basis of colour-magnitude diagram of HDF early-type galax- 
ies. Franceschini et al. (1998) compared model spectra with 
a seven-colour broad band photometry of HDF early-type 
galaxies, and concluded that the major episodes of star for- 
mation building up typical M* galaxies took place during 
1 < z < 4 for the deceleration parameter qo — 0.5 (1 < z < 3 
for qo = 0.15). On the other hand, Bernardi et al. (1998) 
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showed, on the basis of the MgII-cr relation, that nearby 
early-type galaxies in clusters and in the field environment 
consist of the same stellar population, and claimed that the 
bulk of stars in the early-type galaxies had to form at z <; 3 
in the field environment as well as in clusters. 

In this paper, we first examine whether the single burst 
model reproduces both the observed number counts in the 
B, I and K bands and redshift distributions of two samples 
of HDF galaxies selected in the / and K bands. We will find 
that the single-burst model overpredicts galaxies at z 1.4 
in redshift distribution even if the effect of dust extinction 
is taken into account during star formation. Following this 
result, we then introduce models which have longer periods 
of star formation than the single-burst model so that dust 
extinction work longer. We compare these models with the 
observations, and find that with a certain combination of 
parameters there are models whose predictions are consis- 
tent with the observations. Our study can be regarded as 
an extension of the work by Franceschini et al. (1998) in 
which they compared the single burst model and a model 
having a longer star formation period with K-band number 
counts and redshift distribution. We limit our discussion to 
so-called pure luminosity evolution models, because there 
is non-controversial evidence for number evolution or mor- 
phology evolution for early- type galaxies. 

The structure of this paper is as follows. In §2, we de- 
scribe the data used for our analysis, and in §3 we show 
how the local luminosity function is normalised. We com- 
pare the observed number counts and redshift distributions 
with models in §4. Our conclusions are given in §5. The cos- 
mological parameters we have chosen to use throughout this 
paper are (fio, Ao, to) = (0.1, 0, 15 Gyr), unless otherwise 
stated, where Qo is the density parameter, Ao is the cosmo- 
logical constant, and to is the present age of the universe 
(for this choice h = 0.59 where H = lOO/i km s _1 Mpc -1 ). 
Note, however, that we have also examined models for two 
other cosmologies, (f2o, Ao, to) = (1.0, 0, 13 Gyr) and (0.1, 
0.9, 15 Gyr), and found that main results do not change. 



2 DESCRIPTION OF THE DATA 

We use B-, I- and K-band number counts and redshift dis- 
tributions of early-type galaxies. The B-band number count 
in the bright magnitude range is obtained using our mosaic 
CCD camera (Kashikawa et al. 1995) attached to the Las 
Campanas 1-m Telescope in 1995 October and to the 4.2-m 
William Herschel Telescope in 1996 April (4.30 deg 2 in total; 
hereafter MCCD data). The I- and K- band number counts 
in the bright magnitude range are taken from Huang, Cowie 
& Luppino (1998). Number counts in the faint magnitude 
range are taken from work based on the Hubble Space Tele- 
scope (HST) observations (Driver, Windhorst & Griffiths 
1995; Driver et al. 1995; Glazebrook et al. 1995; Casertano 
ct al. 1995; Abraham et al. 1996a,b; Odewahn et al. 1996; 
Driver et al. 1998; Franceschini et al. 1998). 

In deriving number counts of early-type galaxies, mor- 
phological classification is needed. MCCD data are classi- 
fied into early-type and late-type galaxies using the central 
concentration and the average surface brightness of galaxies 
(Doi, Fukugita & Okamura 1993). A part of /-band HST 
data (Abraham et al. 1996a,b) is classified using the central 



concentration and the asymmetry of galaxies. Morphologi- 
cal classification of other data is mostly based on eye-ball 
inspection, though some authors used profile fitting or neu- 
ral network as well (Casertano et al. 1995; Odewahn et al. 
1996; Driver et al. 1998; Franceschini et al. 1998). Morpho- 
logically classified K-band counts are available from only 
two sources (Huang et al. 1998; Franceschini et al. 1998). 

We use redshift distributions of the HDF galaxies se- 
lected in the / and K bands. The J-band sample is taken 
from Driver et al. (1998). There are 47 galaxies in 22 < I < 
26. Only 12 out of the 47 galaxies have spectroscopic red- 
shifts, and photometric redshifts (Fernandez-Soto, Lanzetta 
& Yahil 1999) are used for the remaining 35. The K-band 
sample is taken from Franceschini et al. (1998). There are 
34 galaxies in 16.15 < K < 20.15. Out of the 34 galaxies 14 
have spectroscopic redshifts and 20 are assigned photometric 
redshifts (Franceschini et al. 1998). 



3 NORMALISATION OF LOCAL 
LUMINOSITY FUNCTION 

To predict number counts and redshift distributions, the 
local luminosity function is needed. Recently, the local lu- 
minosity function of early-type galaxies has been derived 
using large redshift surveys. We consider two recent sur- 
veys, SSRS2 (Marzke et al. 1998) and 2dF (Colless 1998). 
The values for M* and a are different between them ((a, 
M* - 5 log h) = (-1.00, -19.37) for SSRS2 and (-0.486, 
-19.46) for 2dF). We adopt the values given in SSRS2, be- 
cause the 2dF results are preliminary. However, we confirm 
that the differences in M* and a between SSRS2 and 2dF 
do not significantly affect our results. The measured val- 
ues for normalisation factor </>* are different by about a fac- 
tor of 2 between SSRS2 (0.0086 /i 3 Mpc~ 3 ) and 2dF (0.0044 
fo 3 Mpc~ 3 ). We do not take these cj>* value. Instead, we de- 
termine 4>* so that the predicted counts at 17.5 < B < 20 
match our MCCD data. The resulting value <f>* = 0.0067 
h Mpc~ 3 , which is ~1.5 times lower than the SSRS2 value, 
is adopted. 



4 ANALYSIS 

4.1 Single burst model of early-type galaxies 

We assume pure luminosity evolution for early-type galaxies 
without the evolution of number or morphology. We com- 
pute galaxy evolution using the spectral synthesis code by 
Kodama & Arimoto (1997) (hereafter KA97). The star for- 
mation rate of galaxies if) is expressed as 

f ^ = TiF M s- (* < *«w) m 
\ V - (t > t gw ) 1 ' 

where tsf is the star formation time scale, t gw is the time 
when galactic wind blows, and M gas is the gas mass. Lu- 
minosity evolution of early-type galaxies, which we are con- 
cerned with, is characterized mostly by two parameters, tsf 
and tgw. KA97 found that the observed colour- magnitude 
relation of early-type galaxies in clusters is well reproduced 
if they adopt tsf =0.1 Gyr, t gw < 0.353 Gyr, and the 
galaxy formation redshift zf — 4.5 (qo = 0.5, h = 0.5). We 
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Figure 1. Redshift distribution for the A'-band sample (16.15 < 
K < 20.15) (histograms). Predictions of the single-burst models 
with zp = 3, 5, 8 and oo are shown by lines. Panels (a) and (b) 
are for no extinction and for the strongest extinction we consider, 
respectively. 



limit our discussion to models with zf > 3, because early- 
type galaxies in clusters are suggested to have zf ;> 3 (e.g., 
Ellis et al. 1997). Hereafter, we call such models with (tsf, 
t gw , zf) = (0.1 Gyr, 0.353 Gyr, zf > 3) the single-burst 
models. For the initial mass function, we assume a power- 
law mass spectrum with a slope of x = 1.10 in the range of 
O.IMq < M < 60A/q. We examine whether the single-burst 
models can reproduce both observed redshift distributions 
and number counts. Figjl|(a) compares the observed red- 
shift distribution of the /f-band sample (solid histogram) 
with predictions by a series of single-burst models with dif- 
ferent zf- Four lines correspond to zf — 3, 5, 8 and oo. It is 
found that the single-burst models overpredict the number 
of galaxies at z £ 1.4 irrespective of zf- In particular, any 
single-burst model predicts a sharp peak near zf, because 
galaxies just after their formation epoch are very bright. 
Such a peak would disappear if we assume that the forma- 
tion redshift of early-type galaxies is not a single value but is 
distributed uniformly in time between zf — 3 and oo. How- 
ever, the overprediction at z £ 1.4 still remains, and these 
models are not consistent with the observations of redshift 
distribution. We also compare predictions with the observed 
redshift distribution of /-band sample, and find similar over- 
predictions. The models shown Fig.[j](a) overpredict the K- 
band counts, but they are marginally consistent with the 
B-band and /-band counts. However, they are ruled out by 
the redshift distributions anyway. 



4.2 Models with interstellar dust 

In order to suppress the overprediction found for the single- 
burst model, we introduce the effect of dust extinction. The 
paucity of galaxies at z Tz 1.4 may not be due to dust ex- 
tinction but to some observational selection effects. The dim- 
ming of surface brightness by a factor of (1 + z) 4 , for exam- 
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Figure 2. The filled circles indicate a combination of parameters 
with which the model reproduces observed number counts and 
redshift distributions for appropriate choice of zp and E(B — V) p , 
while open circles show the parameters which do not account for 
the observations. 



3 15 r(a) 22</ 814 <26 

02 





15 


z F =3 




--- z F =5 




z F =8 




Z = oo 



Figure 3. Redshift distributions (histogram) of the /-band sam- 
ple (a) and the /C-band sample (b). Predictions of models using 
rgp = 0.5 Gyr, t gw = 3.0 Gyr and E(B — V) p = 2 are shown by 
lines. Four lines correspond to zp = 3, 5, 8 and oo. 



pie, would decrease the detection rate of early-type galax- 
ies at high redshifts. To investigate this effect, we examine 
the evolution of surface brightness and effective radius of 
the normal early-type galaxies seen in the local universe on 
the basis of the single-burst models, and find that most of 
these galaxies up to z ~ zf survive the selection criteria 
for surface brightness and radius set for Franceschini et al.'s 
(1998) sample, if they are brighter than K = 20.15 mag, 
the limiting magnitude of Franceschini et al.'s (1998) sam- 
ple. This result is also supported by a similar test using the 
early- type galaxies of Franceschini et al.'s (1998) sample for 
which surface brightness, effective radius and redshift are 
given. Accordingly, the surface brightness selection is not a 
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Figure 4. Differential number counts as a function of apparent magnitude. Panels (a), (b) and (c) are for B, I and K band, respectively. 
Predictions of models using tsf = 0.5 Gyr, t gw = 3.0 Gyr and E(B — V) p = 2 are shown by lines. Four lines correspond to zp = 3, 5, 
8 and oo. Dash-dotted lines are for no evolution models. See the text (§2) for data sources. MCCD data are used to normalise the local 
luminosity function. We do not calculate the counts beyond B > 26 mag because of the lack of reliable /("-correction for redshift above 
z = 2. 



problem. Another possible selection effect comes from type 
classification. Early-type galaxies undergoing intensive star 
formation, where a significant amount of gas exists, might 
be classified as late-type galaxies. Such selection effects need 
further investigation. In this paper, we assume that such se- 
lection effects are negligible. 

Following Vansevicius, Arimoto & Kodaira (1997), we 
express the B-band optical depth of a galaxy as 

r(B) = C/ g Z g (2) 

where f g is the gas fraction of the galaxy, Z s is the metal- 
licity of the galaxy, and C is a constant. The internal ex- 
tinction in the B band is Ab = — 2.5 log[exp{— t(B)}] mags. 
We adopt the extinction curve given in Cardelli, Clayton 
& Mathis (1989). Then, the colour excess is expressed as 
E(B - V) ~ 0.245 t(B). Metallicity Z s increases with time, 
while f g is a decreasing function of time, with f g = at 
t > tgw. Therefore, t(B) (and thus E(B — V)) has a peak 
value. We treat the peak value, E(B — V) p , as a free pa- 
rameter and change it from E(B — V) p = to 2. It should 
be noted, however, that the nominal peak value obtained 
from eq. (13) is not necessarily reached at t < t gvi . Actually, 
E(B — V) = for t > t gw , of course. However, it is the nom- 
inal peak value that we change. E(B — V) p = 2 corresponds 
to extraordinary strong extinction. 

Fig. ^(b) compares the observed redshift distribution of 
Jf-band sample (solid histogram) with predictions by mod- 
els for E(B — V)p = 2 with different zf. It is found that 
there still remains the overprediction. This is because the 
extinction, though very strong, works only in a very short 
period (<, 0.3 Gyr) due to the small values for tsf and i gw . 
The models shown Fig.|l](b) are more or less consistent with 
the number counts. However, they are also ruled out by the 
redshift distributions. 

If larger values are adopted for tsf and t gw , the du- 
ration when extinction works becomes longer. In order to 
investigate this effect, we change tsf and t gw as free param- 
eters in 0.1 < tsf < 5.0 Gyr and t gw > 0.5 Gyr, respec- 



tively, and examine whether there exists a combination of 
tsf and t gw that does not overpredict galaxies at z £ 1.4 
and is consistent with observed number counts. We regard 
a model to be consistent with the observations if the model 
reproduces both observed number counts and redshift dis- 
tributions for appropriate choices of z-p and E(B — V) p . Re- 
sults are shown in Fig.^, where the filled circles indicate 
a combination of parameters which is consistent with the 
observations while open circles show the parameters which 
do not account for the observations. We confirm that the 
combinations indicated by filled circles give B — V colours 
at z = which are consistent with that of local early- 
type galaxies (0.8 < B — V < 1.1; Fukugita, Shimasaku 
& Ichikawa 1995). As an example, we present the predicted 
redshift-distributions using tsf = 0.5 Gyr, t gw = 3.0 Gyr 
and E(B - V)„ = 2 in Fig||. We regard the cases when 
zf — 3 in Fig.H as being consistent with the observation. 
This combination is also consistent with the observed num- 
ber counts as shown in Fig.^. In Fig.^|(a), there are several 
galaxies at z ^ 1.5 in the I band sample, while the models 
predict no galaxies there. We do not think, however, that 
this discrepancy is significant. It is plausible that the amount 
of dust extinction of actual galaxies has various values, un- 
like our models that assume the same dust extinction for 
all galaxies. A small fraction of weaker extinction galaxies 
could explain those galaxies at z ^ 1.5. 

We find from Fig.| that if t S f > 0.5 Gyr and t gw > 1.0 
Gyr are taken, the predicted number counts and redshift 
distributions are consistent with the observations. However, 
E(B — V) p > 1 is needed for all cases. This value is much 
larger than the colour excess of Lyman break galaxies at 
z ~ 3, E(B - V) ~ 0.3, found by Sawicki & Yee (1998). 
In Fig.^, we show the result only for the case of £ gw < 5.0 
Gyr. In fact, we investigate the case that galactic wind does 
not blow, and find that the predicted number counts and 
redshift distributions are consistent with the observations 
as well. However, colours of these galaxies at z = are 
inconsistent with observations. 
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It is interesting to examine whether or not the models 
indicated as the filled circles in Fig^| are consistent with the 
observations of early-type galaxies in clusters as well. We 
find that these models are marginally consistent with the 
observed evolution of luminosity and mean colours of cluster 
early-type galaxies up to z = 1.27, the redshift of the most 
distant cluster at present. (We cannot discuss the slope and 
scatter of the colour-magnitude relations of cluster early- 
type galaxies, because we assume here that all early-type 
galaxies have the same age and metallicity.) At z > 1.27, 
however, these models predict quite different evolution of 
luminosity and colours from that based on the single-burst 
models, due to the effect of dust extinction. 

As seen in Fig.^(c), A"-band number counts of Huang et 
al. (1998) at 17 <; K <, 19 are lower than those of Frances- 
chini et al. (1998) by more than a factor of three. The reason 
for this discrepancy is not clear for us. However, if-band 
number counts of Huang et al. (1998) in bright magnitude, 
K ^ 14, where the effect of galaxy evolution is negligible, are 
also lower than the prediction, though the predictions for B 
and / counts at bright magnitudes are in good agreement 
among different observations. Because we cannot reproduce 
Huang et al.'s (1998) data for if-band counts, we ignore 
them in comparing with models. 



5 CONCLUSION AND DISCUSSION 

We investigate pure luminosity evolution models for early- 
type galaxies which do not exhibit either number den- 
sity change or morphology transition, and examine whether 
these models are consistent with observed number counts 
and redshift distributions. We summarize our findings as 
follows: 

(i) We find that single-burst models are inconsistent with 
observed redshift distributions irrespective of zf (if zf > 3), 
due to the overpredictions of galaxies at z ^ 1.4, even if dust 
extinction which is in effect during star formation is taken 
into account. 

(ii) In order for dust extinction to be more effective, we 
change tsf and t sv , as free parameters over the ranges of 
0.1 < tsf < 5.0 Gyr and t gw > 0.5 Gyr, respectively, and 
examine whether there exists a combination of tsf and r, gw 
that does not overpredict galaxies at z <; 1.4 and is also con- 
sistent with observed number counts. We find that models 
with tsf > 0.5 Gyr and t gw > 1-0 Gyr can be made con- 
sistent with the observed redshift distributions and number 
counts, if we introduce strong extinction (E(B — V) p > 1). 

Our finding that the single burst models cannot express 
the evolution of field early-type galaxies well is not new. 
Franceschini et al. (1998) analysed jf-band number counts 
and redshift distribution for field early-type galaxies, and 
found that models which assume the duration of star for- 
mation to be about 3 Gyr with dust extinction during star 
formation are preferred. Our study not only confirms their 
results but also extends their study: using B-, I- and if -band 
number counts and I- and A"-band redshift distributions, we 
systematically study what combinations of tsf and t gw can 
reproduce the observation. He & Zhang (1999) presented a 
redshift distribution for galaxies limited to 22.5 < bj < 24.0 
and B — K > 5.5, which they considered to be early-type 



galaxies, and compared it with the predictions of single- 
burst models and models with tsf = 1 Gyr. They found 
that these models overpredict the observed redshift distri- 
bution at 2 > 0.8, and suggested that number evolution 
may be essential for early-type galaxies. However, we con- 
firm that using our models with strong dust extinction, the 
overpredictions of the redshift distribution can be removed. 
Driver et al. (1998) analysed J-band number counts and red- 
shift distribution for early-type galaxies and found that the 
single-burst models overpredict galaxies in redshift distribu- 
tion, which is consistent with our results. Im et al. (1999) 
analysed /-band redshift distribution for early-type galax- 
ies and found that observations are consistent with those 
expected from passive luminosity evolution or are only in 
slight disagreement with the nonevolving model. However, 
their results are not in conflict with ours, because the red- 
shift distributions they analysed are at / < 21, where most 
galaxies are at z < 1.0. 

It is true that the pure luminosity evolution models with 
dust extinction are not the only model which reproduces the 
observed number counts and redshift distributions. Mod- 
els with number density change or morphology transition, 
such as those based on hierarchical clustering scenarios (e.g., 
Kauffmann, White & Guiderdoni 1993; Baugh, Cole & Frenk 
1996) may account for the observations. However, we can at 
least conclude that early-type galaxies in the field environ- 
ment have a different evolution history from the single-burst 
models, if the observational selection effects are not signifi- 
cant. 
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